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Along with an apparent convergence in estimates of impact-melt 
volumes [1-3] produced during planetary impact events, intensive 
efforts at deriving scaling relationships for crater dimensions have 
also been yielding results [4], It is now possible to examine a variety 
of phenomena associated with impact-melt production during large 
cratering events and apply them to planetary problems. This contri- 
bution describes a method of combining calculations of impact-melt 
production with crater scaling to investigate the relationships be- 
tween the two. 

Calculations of Impact-Melt V olumes: This study uses melt- 

volume calculations that treat vertical impacts into a flat target; the 
projectile and target are described by a modified, material-specific 
Mumaghan equation of state, the details of which have been 
described elsewhere [3,5]. It does not use the “constant-energy 
shell” assumption of Charters and Summers [6] and Gault and 
Heitowit [7], and yields a closer approximation to the results of 
more complex, finite-difference models [1.2]. In an attempt to 
approximate the off-axis stress decay evident in more detailed 
models of the impact process [1,2,8], the particle velocity in the 
target behind the shock front is assumed to vary as co6* 8, where (1 
is the ratio of target to projectile compression, and 8 is the angle 
from the axis of penetration (equivalent to the direction of the 
velocity vector) to the point of interest, measured at the center of 
flow in the target. Energy is added to the target until the projectile 
decompresses, whereupon the detached shock front is treated as a 
thin region whose intensity decays due to entropy production and 
geometric effects. The results of this procedure in calculating melt 
volumes for the impact of chondrite (simulated by an artificially 
dense basalt) into Hardhat granite and Tahawus anorthosite are 
compared in Fig. 1 to those for anorthosite into anorthosite, using a 
finite -difference model [1]. 
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TABLE 1. Impact-melt volumes at terrestrial craters. 


Crater 

Da 

(km) 

Melt Volume 
(fan*) 

Brent 

3.8 

2x10* 

Zapadnaya 

4.5 

1x10' 

H'ncu 

8 

7x10' 

Kaluga 

15 

8 

Logoiak 

17 

7x10** 

Lappajlivi 

17 

8 

Riel 

24 

2x10'* 

Boltyih 

25 

11 

Mistastin 

28 

20 

W. Clearwater 

32 

80 

Kara 

65 

480 

Manicouagan 

100 

1200 

Popigai 

100 

1750 

Sudbury 

200 

8000* 


‘Values are low probably due to mixed nature of target and presence of 
volatiles in sedimentary strata [9]. 

I Assumes Sudbury Igneous Complex and part of overlying Onaping Forma- 
tion are a coherent melt sheet (10]. 


Melt Volumes and Cavity Scaling: Even though it is far from 

complete, the terrestrial dataset is the best one available for the 
determination of impact-melt volumes associated with large cra- 
ters. Values for these volumes for impacts into crystalline rock have 
been estimated from information culled from the literature (Table 1), 
and are presented in Fig. 2 against transient-cavity diameters , which 
were derived from the observed crater diameters with Croft's [11] 
modification scaling relationship. In order to extrapolate from the 
terrestrial data to planetary events, provision must be made for the 
effects of variations in gravity, impact velocity, and target and 
projectile properties. Calculations were performed for the impact of 
chondritic projectiles into Hardhat granite, covering a velocity 
range of 10 to 50 km s '; chondrites are considered to be represen- 
tative of silicate impactors, and the velocity range falls within that 
allowable for the Earth. The scaling relationship as given by 



Impact Velocity (km s' 1 ) 


Fig. 1. Volumes of pure melt and pure vapor produced for impacts 
of chondrites into granite and anorthosite and for anorthosite into 
anorthosite [ 1 ] . The chondrite calculations were performed with the 
model described above. 



Fig. 2. Comparison between model calculations and actual 
terrestrial data for craters in crystalline rock. Note the displacement 
between the three curves and the observed values. Because they 
were formed in mixed targets, the Rics and Logoisk structures were 
not included in the analyses. 
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Schmidt and Housen [4] was then used to estimate transient-cavity 
diameters for a range of projectile sizes. Since the melt volumes are 
originally expressed in terms of impactor volume, it is a simple 
matter to relate the volume of impact melt to the dimensions of the 
transient cavity as given by the scaling relationship [4], Three such 
curves are also included in Fig. 2. It is immediately apparent that the 
theoretical curves overestimate the quantity of melt by, in some 
cases, more than an order of magnitude. We believe that this reflects 
an underestimate of the transient-cavity dimensions rather than an 
overestimate of melt production for the following reasons: (1) The 
scaling relationship used for cavity dimensions was formulated 
partly on the basis of the final dimensions of craters formed in sand, 
which almost certainly represent adjusted transient cavities. (2) The 
melt-volume estimates are accurate to within a factor of 2, with 
underestimates being equally likely as overestimates, and thus 
cannot account for the differences. (3) Melt ejection could account 
for removal of up to 30% of the total produced at the smallest craters, 
bin will have a vanishingly small effect in the cases of the largest 
craters in the figure. (4) The melt volumes calculated here, as 
evidenced by Fig. 1, are in good agreement with those determined 
from the more complex models. Lacking a detailed physical basis 
for changing the scaling relationship — which, it must be empha- 
sized, combines with the melt calculations to yield a slope that is 
statistically indistinguishable from that of the terrestrial data — the 
model curve is brought into agreement with the terrestrial data 
simply by multiplying the model relationship for 23 km s-' (the rms 
terrestrial impact velocity [12]) by a constant. The resulting rela- 
tionship is 

D tc = 1.39^j 3 D° 7 *V i °- 43 g- <> - 22 (1) 

where D tt is the transient-cavity diameters, p p and p, are the 
projectile and target densities respectively, Vj is the impact velocity, 
and g is the gravitational acceleration, all in cgs units. 

Final Crater Dimensions: Equation (1) can be written for final 
crater dimensions by direct incorporation of Croft's [11] modifica- 
tion scaling relationship, which can be written as 

Dr " DQ ) ' ,8 D{j I * (2) 

in which D R is the final rim-crest diameter and Dq is the diameter 
of the simple-lo -complex transition for the planet (and terrain) in 
question. Substitution for D* into equation (1) and solving for D R 
yields an equation of the form 


Dr = lc^E-J^Dp' 92 Vj 0-51 (3) 

where k is a constant related to g and the value of Dq. Values for Dq 
and k are given in Table 2. Relationships described by equations (1) 
and (3) are used elsewhere in this volume in relating model melt 
volumes to observed characteristics of the terrestrial and planetary 
impact record, and in deriving certain implications of those relation- 
ships for the cratering record. 
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TABLE 2. PbneMpedfic constant! for use in equation (3). 


Planet 

Dq (km) 

k (cm** s"') 

Mercury 

10.0 

2.61x1a 1 

Venus 

4.0 (assumed) 

2J4xl<H 

Earth 

4.0 

2.24x10' 

Moon 

10.9 

3.38x10' 

Mart 

3.1 

3.39x10-' 
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tionships for crater formatio 


Differences in scaling relationships for crater formation and the 
generation of impact melt should lead to a variety of observable 
features and phenomena. These relationships infer that the volume 
of the transient cavity (and final crater) relative to the volume of 
impact melt (and the depth to which melting occurs) decreases as the 
effects of gravity and impact velocity increase. Since planetary 
gravity and impact velocity (Table 1 ) are variables in the calculation 
of cavity and impact-melt volumes [1], the implications of the 
model calculation will vary between planetary bodies; this contri- 
bution will address some of those differences. Details of the model 
calculations of impact-melt generation as a function of impact and 
target physical conditions have been provided elsewhere [1], as 
have attempts to validate the model through ground-truth data on 
melt volumes, shock attenuation, and morphology from terrestrial 
impact craters [2,3], 

Melt Volumes: The volume of impact melt as a function of 

rim-crest diameter is shown in Fig. 1 for typical impact velocities 
at the five terrestrial planets [4] (Table 1). In the calculation of rim- 
crest diameter, a modified version of Schmidt and Housen [5] 
scaling was used to calculate transient-cavity diameters [3], which 
were converted to final rim-crest diameters using the “modification 
scaling” relation of Croft [6] . Chondritic projectiles were used in all 
calculations, and assumed target materials varied by planet (Table 
1). Figure 1 indicates that relative melt volumes at craters of a given 


TABLE 1. Variables used in the calculations of impact melting and 
crater dimensions. All targets were assumed to have a temperature 
of 273 K except for Venus, for which 700 K was used. 


Planet 

Target 

Vj 

(km s ') 

Gravitational 
Acceleration 
(cm s 1 ) 

Mercury 

Anorthosite 

23.6 

370 

Venus 

Diabase 

19.3 

891 

Earth 

Granite 

17.8 

981 

Moon 

Anorthosite 

14.1 

162 

Mars 

Anorthosite 

12.4 

371 
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